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We investigate the ground state structure of the Heisenberg model on the 1/5-depleted square lattice 
for arbitrary values of the first- and second-neighbor exchange couplings. By using a mean- field 
Schwinger-boson approach we present a unified description of the rich ground-state phase diagram, 
which includes the plaquette and dimer resonant-valence-bond phases, long- and short-range Neel 
orders, an incommensurate phase and other magnetic orders with complex magnetic unit cells. We 
also discuss some implications of our results for the experimental realization of this model in the 
CaV40g compound. 



The CaV40g is the first example of a quasi two di- 
mensional magnetic system with a spin gap ill As such, it 
has received much attention recently, with a large part of 
this activity devoted la understanding the mechanism of 
the gap formation.la EJ Most of these works have consid- 
ered the Heisenberg model on the so called CAVO lattice, 
i.e., the 1/5-depleted square lattice of Fig. 1, with both 
first- and second-neighbor interactions (Following the no- 
tation in Ref. 3, we will call J\ and J[ to the nearest- 
neighbor intra and interplaquette couplings respectively, 
and, accordingly, J 2 and J 2 for the second-neighbor in- 
teractions). 
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FIG. 1. The 1/5-depleted square lattice, also called a 
CAVO lattice. We indicate also the first- and second-neighbor 
exchange couplings between vanadium atoms in the CaV^Og 
compound. 



A uariety of methods, including quantum Monte 
Carlo,Bu perturbative and high^temperature series 
expansions,™ exact_diagonalization,H different mean-field 
approximations,Uu, tenormalization group methodslj, 
and other techniquegl3'tll, have been applied to this 
model, leading to the conviction that it has a spin gap for 
J\ <L J\ and/or large enough frustrating second-neighbor 
interactions. Based on simple considerations, most au- 
thors have assumed J{ « J\, « J 2 and J 2 ~ 0.5Ji, 
and, to lowest order, identified the gap as the singlet- 
triplet gap above the plaquette resonant-valence bond 
(PRVB) state obtained from the four nearest-neighbor 
V atoms coupled by J\U However, by using LDA cal- 



culations PicketfEl concluded that the second-neighbor 
interactions might be dominant, with J 2 w 2J\. In 
this case the singlet ground state would correspond to 
weakly-coupled metaplaquettes formed by V atoms con- 
nected by the J 2 bonds. On the other hand, from the 
analysis of neutron inelastic scattering data Kodama et 
aZM estimated J 2 0.1J 2 , J[ ~ J\ ~ O.4J2 in order 
to reproduce the dispersion of the lowest triplet excita- 
tion. Moreover, other authors have suggested that the 
exchange couplings could bejtemperature dependent due 
to spin-phonon interactions n 

All the above mentioned works assumed that a single 
non-degeneratp-jorbital on V atoms is occupied. Marini 
and Khomskiili-3 argued, however, that orbital ordering 
and associated structural distortions determine the mag- 
netic properties of the vanadates CaV n 02 n +i, and pro- 
posed a structure of weakly-coupled dimers along J[ 
bonds for the ground state of CaV40g. Recent work by 
Gelfand and SinghE3 concluded that this proposal is not 
consistent with experimental observations, and that the 
PRVB and metaplaquettes scenarios are both able to pro- 
duce a reasonable explanation for the observed-gap and 
uniform susceptibility. In a very recent workjlj Katoh 
and Imada pointed out serious limitations of the simple 
single-orbital model to reproduce results from neutron 
inelastic scattering. They partially solved these discrep- 
ancies between theoretical predictions and experimental 
observations by considering the effects of orbital degen- 
eracy and orbital order, which lead to an effective spin 
Hamiltonian with exchange couplings strongly dependent 
on patterns of orbital occupancy. 

The above brief summary of previous work on CaV-jOg 
points to the difficulties found to determine the real val- 
ues of the competing couplings involved. Consequently, it 
seems worth to us the investigation of the general ground- 
state structure of the quantum Heisenberg antiferromag- 
net on the CAVO lattice for arbitrary values of J\ , J[ and 
J 2 (we take J 2 = J 2 for simplicity), in order to identify 
possible phases with the properties observed experimen- 
tally. We report here the results of such an investigation, 
which is based on the use of the mean-field Schwinger 
boson approximation. 
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FIG. 2. Ground-state structure of the classical Heisenberg 
model on the CAVO lattice. Full lines indicate the stability 
regions of phases I, Neel and Q = 0. In the shaded area 
an incommensurate phase is stabilized in the thermodynamic 
limit. The dashed lines indicates the stable phases I, II, Neel 
and Q = for a cluster of 32 spins. 



Let us start by considering first the ground state of 
the classical Heisenberg model on the lattice of Fig. 1. 
A standard analysis^ shows that even in this case the 
model turns out to be very rich, as shown in Fig. 2. 
The stability regions of the different phases are indicated 
by full lines; we also show, for further use, the stable 
phases in lattices of up to 32 spins (dashed lines). The 
Q = and Neel phases have antiferromagnetic order 
in the plaquettes which repeats itself along the trans- 
lation vectors S± = (2, l)a, 62 = (— l,2)a with magnetic 
wavevectors Q = (0, 0) and Q = (n, n) respectively (a 
is the nearest neighbor V-V distance). Phases I and II 
have both a 8-spin magnetic unit cell, translation vectors 
Si = (l,3)a, 62 — (— 3, l)a and a magnetic wavevector 
Q = (7r,7r), with the corresponding spin arrangements 
in the unit cells shown in Fig. 3a, b. The angle if veri- 
fies tan^ = — J{/2J2 in phase I and tamp = —2J\/J[ in 
phase II. Finally, the shaded area is the stability region of 
an incommensurate phase very difficult to characterize. 
This makes the consideration of the quantum nature of 
the spins quite involved in this phase (see below). 
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FIG. 3. Magnetic unit cells and displacement vectors <5i 
and 82 corresponding to a) Phase I and b) Phase II (see main 
text). 



We must now correct the classical diagram of Fig. 2 
by introducing the effects of quantum fluctuations. In or- 
der tO|dp so, we will use the mean-field Schwinger boson 
theoryL3, which has already been applied to the CAVQ 
lattice without frustrating second-neighbor interactions.0 
This approach is suitable for the problem at hand since 
it allows to treat magnetic and non-magnetic quantum 
phases in an unified way. Moreover, despite the fact 
that the Schwinger boson theory is known to overes- 
timate the stability of magnetic phases, it most often 
provides a qualitatively reliable phase structure. On 
the other hand, its predictions can be systematically 
corrected by the inclusion of fluctuations in the mean- 
field order parameters which has been shown to bring 
Schwinger boson results well in line with exact values 
on finite lattices and with the predictions of other more 
accurate methods in the thermodynamic limit. 

Using the standard representation of spin operators in 
terms of two boson fields linked by the local boson num- 
ber restrictions, quantum corrections to the phase dia- 
gram of Fig. 2 can be obtained by simply choosing the 
right initial structure of the ©eder parameters in solving 
the self-consistent equations.Ea Although this is a mat- 
ter of a straightforward energy minimization, the large 
number of spins in the magnetic Brillouin zone makes 
the calculations involved and computationally demand- 
ing. In addition, since we have no simple characteri- 
zation of the incommensurate phase, in this region we 
had to rely on an extensive search for global solutions of 
the minimization equations. The strategy used was first 
to determine the classical ground state in clusters of up 
to 64 spins, and then searching for the absolute energy 
minimum in the quantum incommensurate phase start- 
ing from these large magnetic cells. In particular, up to 
32 spins the stable classical phases are those showed in 
Fig. 2 with dashed lines. In addition to the commensu- 
rate phases above described, there is another phase sta- 
ble only along the dashed line that separates phases I and 
II; it has a 4-spin magnetic cell and magnetic wavevector 
Q = (0,7r). This classical phase is particularly interest- 
ing as a starting point for the search of order parameters 
in the quantum incommensurate phase. This is so be- 
cause no pair of classical vectors in the magnetic motif 
are collinear or anticollinear, which means that the fer- 
romagnetic and antiferromagnetic channels are both ac- 
tive in the minimization. Furthermore, for clusters with 
more than 32 spins the resulting classical stable phases 
are incommensurate (span the cluster). We stress that 
this investigation of the incommensurate phase is allowed 
only because of the particular properties of the Schwinger 
boson theory, since other semiclassical approaches (spin 
wave theory, for instance) require the precise knowledge 
of the classical ground-state structure. Moreover, the in- 
clusion of quantum (non-magnetic) phases in the diagram 
of Fig. 2 is possible because the Schwinger bosons pre- 
serve the rotational invariance of the Heisenberg Hamil- 
tonian when there is no condensation (magnetization) in 
the system. 
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FIG. 4. Ground-state structure of the quantum Heisenberg 
model on the CAVO lattice. Full thick(thin) lines indicate 
first (second)-order phase transitions; the dashed line sepa- 
rates long-range and short-range Neel orders. See main text 
for a description of the different phases. 



After the consideration of quantum fluctuations at 
mean-field order, the ground state has the complex struc- 
ture shown in Fig. 4. The Q = 0-PRVB transition is 
weak first order, while the incommensurate-short-range 
Neel transition might be a very weak one instead of the 
second-order transition indicated. The Chain phase cor- 
responds to disconnected zig-zag chains along J{and Ji 
bonds, with short-range antiferromagnetic correlations 
between nearest-neighbor spins in the chains (this phase 
is essentially the ground state proposed by Marini and 
KhomskiO) . The PRVB and dimer phases are the stan- 
dard singlet phases discussed in previous works, con- 
sisting respectively of isolated plaquettes made up by 
Ji bonds and the J[ dimer covering. For J' 2 3> J[ 
phase I gees to Pickett's weakly-coupled metaplaquette 
scenario.E3 Notice that in this limit the angle tp between 
metaplaquettes becomes 7r/2, instead of or tt expected 
from order-from-disorder ideas. Of course, this phase di- 
agram might be strongly modified by order-parameter 
fluctuations, but it gives an idea of the model complex- 
ity and becomes a starting point for these much more 
involved calculations. 

Notice that the most studied region in connection with 
the CaV40g compound (J{ rj Ji, J2 ~ 0.5Ji) corre- 
sponds to a region where several distinct phases, merge. 
On the other hand, Pickett's proposal J2 « 2 J1E2 is well 
inside the phase I. Then, it is of interest to discuss in 
more detail what happens along the line J{ = J\ and 
arbitrary J2. In Fig. 5 we see the behavior of the energy, 
with the continuous transition from the short-range Neel 
order to the incommensurate phase at Ji — 0.55Ji, and 
then a first-order transition from this last phase to phase 
I at J2 — 0.62 Ji. Fig. 6 shows the corresponding magne- 
tization curve, which indicates a robust magnetic order in 
phase I and a rather peculiar behavior in the incommen- 



surate phase. In this last phase the magnetization m has 
an abrupt drop to zero near the transition to the short- 
range Neel order, which is in agreement with the smooth 
behavior shown by the energy However, the incommen- 
surate magnetic vector seems to have a discontinuity be- 
fore reaching (tt, tt) (see Fig. 7). Nevertheless, this could 
be only a problem with the numerical accuracy, since 
the convergence of the selfconsistent equations becomes 
more and more difficult near the transition. Another re- 
markable characteristic of the magnetization curve is the 
local minimum in the incommensurate phase. Since the 
behavior of Q is smooth across these points, this could 
indicate that the spins inside the magnetic unit cell rear- 
range without changing the incommensurate wavevector. 
Such an explanation is supported by the small discon- 
tinuities we observed in the short-range correlations on 
both sides near the deep local minimum of m. Finally, 
in Fig. 6 the dashed and full lines in the incommensu- 
rate phase give the magnetization of different sublattices. 
These two sublattices correspond to the zigzagging paths 
along J{ and J2 bonds, in one direction and in the direc- 
tion perpendicular to it respectively. However, in this 
case the paths are not disconnected like in the Chain 
phase, since correlations in J\ bonds do not vanish. In 
passing, we note that to obtain this solution from the self- 
consistent equations one must allow for different values of 
the Lagrange multipliers that impose the boson-number 
restriction on inequivalcnt sites of the magnetic unit cell. 
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FIG. 5. Ground-state energy along the line J[ = Ji in pa- 
rameter space. The arrows at J2/J1 — 0.55 and 0.62 indicate 
the points where there are phase transitions. 
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FIG. 6. Magnetization m as a function of J2/J1 for 
J[ = Ji. The full and dashed lines in the incommensurate 
phase correspond to the magnetization of two different sub- 
lattices. 



instance, in the J\ — J2 model on the square lattice&3). 
In the absence of frustration, the wavenumber of the low- 
est energy excitation (corresponding to the singlet-triplet 
gap in the PRVB scenario) has been predicted from-Ber- 
turbation expansionsB and variational Monte Carldld to 
be (7r,7r), and it shifts to incommensurate wavenumbers 
when J2 becomes relatively large. Our results for the 
Neel, PRVB and I phases indicate also a (ir, ir) wavec- 
tor for the lowest excitations, while the incommensurate 
phase has a Q that changes with J2 as shown in Fig. 
7. Howawr, recent neutron inelastic scattering data for 
CaV40gll3 are consistent with Q = (0, 0), like in a disor- 
dered Q = phase. Finally, we are currently exploring 
the possibility of introducing fluctuations in the order pa- 
rameters, although these are quite involved calculations. 
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FIG. 7. The magnetic wavevector Q/n in the incommen- 
surate phase as a function of J2/J1 for J[ = J\. 



In conclusion, we have determined the general ground- 
state phase structure of the frustrated Heisenberg model 
on the CAVO lattice by means of the Schwinger boson 
mean-field theory. In particular, near the most studied 
parameter region for the CaX-^Og compound we iden- 
tified several competing phases. Although the relative 
stability of these phases can be strongly modified by the 
introduction of Gaussian fluctuations in the order param- 
eters, it is nonetheless quite satisfactory to have a unified 
picture of the phase diagram at mean-field order. More- 
over, it provides a necessary starting point for the more 
involved one-loop calculations. On the other hand, de- 
spite the fact that some of the phases have long-range or- 
der -and, consequently, no gap-, the reduced magnetiza- 
tion values suggest that fluctuation corrections might eas- 
ily destroy the magnetic order and produce the gapped 
phase observed experimentally (like what happens, for 
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